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Abstract 


\ 

The  report  examines  the  feasibility  of  neutralizing  the  potential 
difference  that  occurs  on  a  spacecraft  by  emitting  ions  from  the  positively 
charged  surface  region  of  the  spacecraft  and  lefting  this  current  impact 
on  the  negatively  charged  surface  region.  The  investigation  was  performed 
by  making  two  different  simplifying  assumptions  about  the  geometry  of  the 
problem.  In  the  first  case  we  approximated  the  geometry  by  two  conducting 
infinite  half  planes.  In  the  second  case  we  approximated  the  geometry  by 
two  finite  width  constant  charge  density  plates.  In  both  cases  we  investi¬ 
gated  the  dependence  of  the  particle  trajectories  on  the  origination, 
velocity,  and  direction  of  emission  of  the  particles.  Two  generalizations 
emerge  from  the  study:  1)  Shallow  launch  angles  give  more  favorable 
trajectories.  2)  For  a  given  launch  angle  there  is  an  optimum  energy  which 
yields  a  minimum  impact  distance. . 

A - 


iii 


I. 


INTRODUCTION 


The  investigation  completed  for  this  report  studied  the  feasibility  of 
neutralizing  the  potential  difference  which  occurs  on  the  surface  of  a  space¬ 
craft.  It  has  been  proposed  that  neutralization  can  be  accomplished  by  emitting 
ions  from  the  positively  charged  surface  region  of  the  spacecraft  and  by  letting 
this  current  impact  on  the  negatively  charged  surface  region.  In  order  to  in¬ 
vestigate  the  main  features  of  the  proposed  scheme  we  made  some  simplifying 
assumptions  about  the  geometry  of  the  problem.  In  the  first  case  we  approximated 
the  geometry  by  two  conducting  infinite  half  planes  and  investigated  the  trajec¬ 
tories  of  a  positive  ion  emitted  at  some  point  on  the  positively  charged  half 
plane  with  given  initial  velocity  and  given  direction  of  emission. 

In  the  second  case  we  approximated  the  geometry  by  two  constant  charge 
density  finite  width  strips  and  investigated  the  trajectories  of  a  positive 
ion  emitted  at  some  point  on  the  positively  charged  strip  with  respect  to 
origination,  velocity  and  direction.  The  second  case  probably  provides  a  more 
realistic  description  of  actual  satellite  conditions. 

H.  POTENTIAL  DUE  TO  TWO  CONDUCTING  INFINITE  HALF  PLANES 

A.  Statement  of  problem 

This  case  was  already  discussed  in  Scientific  Report  #1.  We  provide  here 
a  review  of  the  problem  along  with  more  complete  solutions  and  conclusions. 

We  consider  the  plane  y«*o  and  assume  a  cut  on  this  plane  along  the  z  axis. 

Let  there  be  a  given  potential  difference  V  between  the  two  half  planes  x  >  o  and 
x  <  o.  The  equations  of  motio.i  of  a  charged  particle  in  the  electric  field 
produced  in  such  a  configuration  are 
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2  2 
x  +y 


2 .  2 
x  +y 


Motion  in  the  z  direction  is  trivial  and  not  related  to  the  problem.  The 
remaining  equations  can  be  cast  in  a  dimensionless  form 


dt2  x2+y2 


d  y  _  x 
dt2  x2+y2 


with  the  initial  conditions  that  at  time  t=o 


x  =  1 


y  =  o 


v  =  v 
y  xo 


and  where  the  initial  velocities  v  and  v  are  in  units  of  A  ■  - 

xo  yo  mil 


B.  Numerical  Solution 

Since  no  analytical  solutions  to  the  above  equations  could  be  found  they  were 

solved  numerically  using  a  modified  version  of  scientific  subroutine  DHPCG  based 

on  Hamming's  Modified  Predictor  Corrector  Method.  Table  1  lists  the  values  of  xA 

2  2~ 

v  +v  (the  total  initial 
i  x  y 

velocity)  and  the  initial  launch  angle  6.  Table  2  lists  the  values  of  Y  (the 

max 

1  2  2~ 

maximum  height  attained  by  the  first  particle  for  various  values  of  v  “\v  +v 

1  x  y 

and  the  initial  launch  angle.  Data  for  the  blank  areas  of  both  tables  were  not 


obtainable  in  a  reasonable  amount  of  computer  iterations  due  to  their  large  size. 
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The  discrepenci.es  between  these  results  and  the  earlier  ones  given  in  the 
previous  report  for  larger  values  of  v  are  due  to  defects  in  the  earlier  computer 
program.  Figure  1  shows  the  dependence  of  x^  (the  impact  distance)  on  vt  (the 
total  initial  velocity)  and  9  (the  launch  angle).  Figure  2  shows  the  dependence 
of  y  (the  maximum  height)  on  v  (the  total  initial  velocity)  and  6 (the  launch 
angle).  The  figures  show  that  both  xfc  and  y^^  have  a  similar  dependence  on  vt 
and  a  Note  also  that  for  a  given  launch  angle  (except:© =0)  both  x^  and  yfflax  have 


/gv 

—  )  • 

c  l  —  mir 

C.  Asymptotic  Result 

As  can  be  seen  in  figure  1  for  0°  launch  angle  the  impact  distance  x^  tends 
to  0  for  large  initial  velocities.  This  can  be  demonstrated  analytically  by 
approximating  the  basic  equations 


7t.  2 
x  +y 


For  large  initial  velocity  v^o  and  0  launch  angle  the  trajectory  height,  x*, 
is  expected  to  remain  close  to  zero.  Approximating  y»0in  the  above  equations 
leads  to 

x  ■  0  which  yields  x  =  1  +  v  t 

xo 

..  1 

y  “  1+v  t 
xo 

This  equation  may  be  integrated  immediately  and  thus  gives  for  y 


y  -  ^  l»<Hv„0t). 

XO 


1 
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The  equation  for  y  may  be  obtained  by  integrating  again 
y  =  — — «  |(l+v  t  )ln(l+v  t)  -  (1+v  t)+l 

^  2  L  xo  x0  xo  -I  « 

xo 

To  find  the  impact  aistance  x*,  we  set  y=o  and  substitute  for  x  *=  1+v  t 

0  =  xlnx-x+l 

or  lnx  »  1  -  — 

x 

To  find  solutions  for  negative  x  we  can  rewrite  the  equation  as 

lnx  =  1  +  — 
x 

and  solve  for  positive  x. 

This  yields  x  =  3.59112  as  a  solution. 

Therefore  we  see  that  the  launch  angle  0=0°  curve  in  figure  1  has  as  an 
asymptote  the  line  x=3.59  in  the  limit  of  large  initial  velocity.  This  impact 
distance  x*  =  3.59  would  be  the  shortest  obtainable  impact  distance  under  the 
most  favorable  conditions  in  this  geometrical  approximation.  The  data  in  Table  1 
verify  this  conclusion. 

III.  POTENTIAL  DUE  TO  FINITE  WIDTH  CONSTANT  CHARGE  DENSITY  PLATES 
A.  Statement  of  the  Problem 

We  consider  the  plane  y=0  and  assume  a  cut  on  this  plane  along  the  Z  axis.  Let 
there  be  a  constant  charge  desnity  in  the  two  infinite  length  strips  -a<x<o  and 
o<x<a.  One  can  calculate  the  potential,  for  this  geometry  by  starting  with  the 
potential  for  a  line  charge  and  integrating  over  the  area  of  the  strips  with  the 


result 
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v  -  -  j  x  In  [  t(^a).2+y^.LLx+al2±z!}]  +  a  In  [  fet §£  +  jl  1 
*  i  L  (x2+y2) 2  J  L  (x-a) 2+y2  J 

+  2y  ^  tan-1  (^>  +  tan"1  (^)  -  2  tan"1  j  | 


The  electric  field  follows  directly  from  the  potential. 


„  V  r  C(x-a) 2-fy2]  [(x*a)24-y2]  1 

11  L  (x2+y2)2  J 

=  £  2  | tan"1  (m)  +  tan'1  (^)-2  tan"1/'  *)\ 

*  \  \y  /  y  \y J) 


E  =  o 
z 


The  equations  of  motion  are  therefore 


mx 


<*V)2 


..  eV  I 
my  =  —  2  \  tan 


mO  =  o 


1  tan"1 


-1  (£±5)  _2  tan“*  (i) 


"(t)  I 


Again  motion  in  the  z  direction  is  trivial  and  not  related  to  the  problem. 
The  remaining  equations  can  be  case  in  a  dimensionless  form 

x  =  ln(l^-al!t|!]i(|tal2^L2l  j 

I  (x2+y2)2  J 


y  -  2  [tan’1  (*=*)  +  tan'1  (**)-  2  tan'1  (f)  j 


with  the  initial  conditions  that  at  time  t“o 
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x  *  x 

o 


V  =  V 
X  xo 


y  =  o 


V  =  V 

y  yo 


where  the  initial  velocities  v  and  v 

xo  yo 


are  in  units  of 


/  eV 
mir  ' 


B.  Numerical  Solution 

These  equations  were  solved  numerically  using  a  modified  version  of  scientific 

subroutine  DHPCG  based  on  Hamming's  Modified  Predictor  Corrector  Method.  (For  the 

sake  of  the  numerical  calculations  we  assumed  the  width  of  the  strips  a=10.)  Table 

3  lists  the  values  of  x^  (the  impact  distance)  for  various  values  of  v^_  =  Jv ^  +  v ^ 

(the  total  initial  velocity)  and  8  (the  launch  angle)  with  the  launch  point  xq  =  -1. 

Table  4  lists  the  values  of  y  (the  maximum  height)  for  various  values  of  V  and  0 

max  t 

again  for  xo=-l.  Figures  3  and  4  show  the  variation  of  the  test  particle  trajectories 
with  different  launch  positions  for  vt**o.  Note  that  above  xq  =  -7-. 25  the  trajec¬ 
tories  no  longer  bend  to  the  right  and  that  even  for  xq= -1  the  impact  is  already 
beyond  the  width  of  the  strips  (a=10).  Figure  5  shows  the  dependence  of  the  impact 
distance,  x.,  on  the  total  initial  velocity  V  and  launch  angle  0.  Figure  6  shows 
the  dependence  of  the  maximum  height  y  ^  on  the  total  initial  velocity  and 
launch  angle  0.  As  in  the  previous  case  both  x^  and  yjnax  have  a  very  similar  de¬ 
pendence  on  these  parameters.  In  this  case  both  x^  and  ymax  have  minimums  for  the 
total  initial  velocity  in  the  range  . 5<Vt<2  for  launch  angles  below  40°  only.  Figures 
7,  8,  and  9  give  a  more  detailed  picture  of  the  dependence  of  the  impact  dis¬ 
tance  and  of  the  maximum  height  of  the  trajectory  on  the  parameters.  Figure  10 
shows  the  dependence  of  the  impact  distance,  x*,  and  maximum  height,  y  ,  on  the 
initial  launch  position,  xq,  for  various  angles  9.  Note  that  a  small  increase  in 
Xq  can  lead  to  vary  large  increases  in  xA  and  ymax»  This  dependence  was  not  pre¬ 
sent  in  the  first  case. 
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IV.  CONCLUSIONS 

We  have  examined  particle  trajectories  in  two  simple  electric  field 
configurations.  In  both  cases  we  found  trajectories  which  struck  the  target. 
Two  generalizations  emerge  from  this  study:  1)  shallower  launch  angles 
lead  to  more  favorable  trajectories  and  2)  for  a  given  launch  angle  there  is 
an  optimum  energy  which  yields  a  minimum  impact  distance.  Whether  these 
generalizations  remain  true  in  more  complicated  geometries  remains  to  be 
investigated. 
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CAPTIONS 

1.  Conducting  plates:  plot  of  impact  distance  xA  as  a  function 
of  the  total  initial  velocity  v  for  various  launch  angles  0. 

2.  Conducting  plates:  plot  of  the  maximum  height  of  the  trajectory 

y  as  a  function  of  the  total  initial  velocity  v  for  various 
max  J  t 

launch  angles  0. 

3.  Dielectric  plates:  plot  of  trajectories  as  a  function  of  launch 
position  xq  (v  =0) 

4.  Dielectric  plates:  plot  of  trajectories  as  a  function  of  launch 
position  xq  (v  =0) 

5.  Dielectric  plates:  plot  of  impact  of  distance  x^  as  a  function 
of  the  total  initial  velocity  v£  for  various  launch  angles  0. 


6.  Dielectric  plates:  plot  of  the  maximum  height  of  the  trajectory 

y  as  a  function  of  the  total  initial  velocity  v  for  various 
max  t 

launch  angles  0  (xq  =  -1) 

7.  Dielectric  plates:  plots  of  x*  and  ymax  as  functions  of  the  initial 
horizontal  velocity  (xq  =1,  0=0) 

8.  Dielectric  plates:  plots  of  x.  and  y  as  functions  of  the  initial 
vertical  velocity  (xq  =  1,  0  =  90) 

9.  Dielectric  plates:  plots  of  x.  and  y  as  functions  of  the  launch 

r  r  *  max 

angle  0  (xq  =  1,  v£  =  1.0) 

10.  Dielectric  plates:  plots  of  x.  and  y  as  function  of  the 

*  max 

initial  launch  position  xq  for  various  angles  8.  (v  =  1.0) 
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TABLE  1 
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